UNCLASSIFIED  TR-0077(2940)-3  SAMSO-TR-77-154  NL 


NATIONAL  BUREAU  OF  STANDARDS 

MICROCOPY  RESOLUTION  TEST  CHART 


i*n  m* 

rvJ  nu. 

BBC  FILE  COPY  AO  A 0 4 3 3 3 9 


Simple  Model  of  a Chain-Reaction 
Pulsed  HF  Laser 


Department  of  Mechanical  Engineering 
Michigan  State  University 
and 

Aerophysics  Laboratory 
The  Ivan  A.  Getting  Laboratories 
The  Aerospace  Corporation 
El  Segundo,  Calif.  90245 


29  July  1977 


Interim  Report 

D D C 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


^ B 

Prepared  for 

SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
AIR  FORCE  SYSTEMS  COMMAND 
Los  Angeles  Air  Force  Station 
P.O.  Box  92960,  Worldway  Postal  Center 
Los  Angeles,  Calif.  90009 


_iiMK 


,K' 


4 


r 


This  report  was  submitted  by  The  Aerospace  Corporation,  El  Segundo, 
^ CA  90245,  under  Contract  F0470  1-76-C-0077  with  the  Space  and  Missile 

Systems  Organization,  Deputy  for  Advanced  Space  Programs,  P.O. 

Box  92960,  Worldway  Postal  Center,  Los  Angeles,  CA  90009.  It  was  re- 
viewed and  approved  for  The  Aerospace  Corporation  by  W.  R.  Warren,  Jr.  , 
Director,  Aerophysics  Laboratory.  Lieutenant  A.  G.  Fernandez,  SAMSO/ 
YAPT,  was  the  project  officer  for  Advanced  Space  Programs. 

This  report  has  been  reviewed  by  the  Information  Office  (OI)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publica- 
tion. Publication  of  this  report  does  not  constitute  Air  Force  approval  of 
the  report's  findings  or  conclusions.  It  is  published  only  for  the  exchange 
and  stimulation  of  ideas. 


FOR  THE  COMMANDER 


Asst.  Deputy  for  Advanced  Space 
Programs 


UNCLASSIFIED 


SECURITY  Cl^A&SlFlCATION  OR  THIS  PAGE  (m\tn  D ata  End'd) 


READ  INSTRUCTIONS 


EPORT  DOCUMENTATION  PAGE 


BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CATALOG  NUMBER 


FgOVT  ACCESSION  NO 


S AMSO-iTR  - 77-  1 54 


5 TYPE  OF  REPORT  ft  PERIOD  COVERED 


*.  TITLE  («nd  SubrltlaJ 


SIMPLE  MODEL  OF  A CHAIN -R EAC TION 
PULSED  HF  LASER 


S.  PERFORMING  0 RG.  R&PflfiT  NUMBER 

TR-OO77(2940)-3 

ir  contract  or  grant  numbers; 


7 AuTMOR(J) 


F04701-76-C-0077 


Ronald  L.  Kerber  and  James  S.  Whittier 


PROGRAM  ELEMENT.  PROJECT 
AREA  ft  WORK  UNIT  NUMBERS 


9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


The  Aerospace  Corporation 
El  Segundo,  Calif.  90245 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Space  and  Missile  Systems  Organization  // 

Air  Force  Systems  Command  — " 

Los  Angeles,  Calif.  90009 

14  MONITORING  AGENCY  NAME  ft  ADDRCSSJY/  different  trom  Controlling  Otllca) 


T3.  NUMBER  OF  PAGES 


IS.  SECURITY  CLASS,  (ol  thla  raport) 


DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  Distribution  statement  (ot  this  Raport) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (ot  tha  abatract  antarad  In  Block  20,  It  dlttaranl  from  Raport) 


KEY  WORDS  (Continue  on  reverae  aide  it  necessary  and  identify  by  blo.k  number) 


Pulsed  Chemical  Laser 
Chemical  Laser  Model 
Chain-Reaction  Laser 


ABSTRACT  (Continue  on  reverae  aide  It  neceeaary  and  Identify  by  block  number) 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWh an  Data  Entarad) 


UNCLASSIFIED 


CONTENTS 


PREFACE  1 

I.  INTRODUCTION  5 

II.  THEORY  9 

A.  Model  Formulation  9 

B.  Approximate  Model  Solution  

C.  Comparison  with  Comprehensive  Model 19 

III.  PARAMETRIC  STUDY  OFLASER  PERFORMANCE  23 


A.  Threshold  23 

B.  Initiation  Level  25 

C.  Diluents  25 

D.  Initial  H^  Concentration  28 

E.  Initial  Temperature  28 

IV.  COMPARISON  WITH  EXPERIMENT  31 

V.  CONCLUSIONS  3 3 

APPENDIX  A 35 

APPENDIX  B 37 

REFERENCES  39 


-3- 


CONTENTS 


PREFACE  1 

I.  INTRODUCTION  5 

II.  THEORY  9 

A.  Model  Formulation  9 

B.  Approximate  Model  Solution  16 

C.  Comparison  with  Comprehensive  Model 19 

III.  PARAMETRIC  STUDY  OFLASER  PERFORMANCE  23 

A.  Threshold  23 

B.  Initiation  Level  25 

C.  Diluents  25 

D.  Initial  Concentration  28 

E.  Initial  Temperature  28 

IV.  COMPARISON  WITH  EXPERIMENT  31 

V.  CONCLUSIONS  33 

APPENDIX  A 35 

APPENDIX  B 37 

REFERENCES  39 


-3- 


i5T'£o.“.DI  (1  P^E  EL-NK-KCT  FILMED 


FIGURES 


1.  Comparison  of  Predicted  Power  Time  Histories  for 

F + H = 0.  1 (a),  0.01  (b),  and  0.  001  (c) 20 

2.  Effect  of  Cavity  Threshold  for  Selected  Levels  of 

Initiation 24 

3.  Effect  of  Level  of  Initiation  for  Selected  Values  of  26 

4.  Effect  of  He  Concentration  for  Selected  Values  of 

Initial  Pressure  27 

5.  Effect  of  H£  Concentration  for  Selected  Levels  of 

Initiation  29 

6.  Effect  of  Initial  Temperature  at  Selected  Values  of 

Initial  Pressure  30 

7.  Comparison  of  Theory  and  Experimental  Results  of 

Hofland  et  al.6  at  800  Torr  32 


TABLES 


I.  Rate  Coefficients  for  H^  + F^  Chemical  Laser 11 

II.  Band  Energy  Level  Spacings 14 

III.  Comparison  of  Pulse  Energy  as  Predicted  by  the 

Simple  and  the  Comprehensive  Model 22 


-4  - 


j. 


.. 


I.  INTRODUCTION 


Interest  in  H 2 + F0  chain- reaction  pulsed  lasers  has  been  increased  by 
steady  experimental  progress  with  electrically  initiated  devices.  Electrical 
efficiencies  greater  than  100%  have  been  achieved,  and  large  output  energies , 
up  to  2360  J,  have  been  achieved  in  electrically  initiated  devices.  These 

developments  have  inspired  increased  interest  in  modeling  of  + F ^ chain 
lasers  for  better  understanding  of  performance  and  for  guidance  in  the  design 
of  improved  devices.  This  report  presents  a simplified  rate  equation  model 
that,  in  comparison  with  earlier  models,  provides  great  economy  of  prediction 
of  laser  power  and  energy  with  little  sacrifice  in  accuracy.  Our  model  has 
some  ad  hoc  features.  However,  we  believe  this  to  be  amply  justified  (1)  be- 
cause the  H_  + F?  laser  is  of  growing  importance,  and  (2)  because  of  the 
computational  economy  that  is  gained. 

Modeling  of  the  pulsed  H7  + F?  laser  was  pioneered  by  Emanuel  and 
7 11  L L 

co-workers.  ' Emanuel  formulated  a "constant- gain " rate  equation  model 

o 

that,  when  solved  numerically,  provides  predictions  of  instantaneous  species 

concentrations,  temperature,  and  the  laser  power  spectrum  of  each  vibrational 

band  as  influenced  by  rates  of  chemical  pumping  and  heating,  molecular  energy 

transfer  and  deactivation,  dissociation  and  recombination,  and  stimulated 

emission.  A study  of  the  sensitivity  of  model  predictions  to  many  parameters 

was  carried  out.  An  early  comparison  of  Emanuel's  model  with  pulsed  laser 

experiments  was  made  by  Suchard  et  al.  , who  found  good  agreement  with 

the  time-to-laser  threshold  and  pulse  shape  of  their  photolytically  initiated 

1 2 

laser.  A related  earlier  study  by  Chester  and  Hess  showed  good  agreement 
between  their  photo-initiated  experiments  and  their  small-signal-gain  rate 
equation  model.  More  recently,  predictions  of  Emanuel's  model  have  been 

13 

compared  with  discharge  laser  experiments;  fair  agreement  was  found  for 
laser  power  and  energy. 


Pulsed  H2  + F2  laser  modeling  research  is  progressing  along  three 
distinct,  but  interacting,  lines.  First,  better  agreement  between  model 
predictions  and  measured  laser  power  and  energy  is  being  sought  through 
kinetic  rate  adjustments,  particularly  those  of  deactivation  processes, 
through  incorporation  of  ionic  species  in  the  laser  kinetics,  13  as  well  as 
through  improved  experiments.  Second,  more  realistic  predictions  of  laser- 
output- spectra  versus  time  as  well  as  modest  refinements  in  energy  predic- 
tions are  being  sought  with  models  that  account  for  cavity  transients  and 
rotational  nonequilibrium  phenomena.  14~  16  The  present  report  belongs  to  a 
third  line  of  investigation  in  which  economical  predictions  are  sought  of 
gross  laser  power  and  energy  through  simplification  of  the  comprehensive 
rate  equation  approach.  It  is  expected  that  the  paper  will  contribute  to  the 
first  two  lines  of  research  as  follows:  The  relative  importance  of  various 
phenomena  is  clarified  in  constructing  the  simple  model  and  in  comparing 
its  predictions  with  those  of  comprehensive  models.  Examination  of  the 
suitability  of  proposed  rate  changes  should  be  facilitated  because  many  para- 
metric combinations  can  be  considered  cheaply.  For  example,  theoretical 
predictions  for  comparison  with  the  entirety  of  available  laser  performance 
data  could  easily  be  repeated  for  a large  number  of  tentative  choices  of 
rate-constant  values.  Moreover,  because  of  its  economy,  the  model  can 
also  be  useful  as  a "building  block"  in  engineering  systems  analysis  studies 

where  laser  size,  weight,  and  performance  are  optimized  for  specific 
applications. 

Previous  work  of  Skifstad  on  a simple  model  for  the  H + F laser 
concentrated  on  the  derivation  of  closed-form  analytical  results;  consequently, 
for  tractability,  his  analysis  assumed  constant  temperature.  17  In  practice, 
such  analytical  results  usually  serve  only  as  a qualitative  guide  to  laser 
performance  as  there  is  generally  a large  temperature  rise  during  the  laser 
pulse  (or  downstream  in  a cw  laser).  In  the  present  work,  we  are  seeking 
accurate  predictions  for  a wide  range  of  conditions;  hence,  we  account  for 
temperature  variation. 


L 


Our  analysis  (Section  II)  stems  from  earlier  work  on  simple  models 
for  the  DF/C02  and  F + pulsed  lasers  18,  ^ In  striving  for  accurate 
prediction  in  the  more  complicated  H2  + F^  pulsed  laser  system,  however, 
we  incorporate  features  that  were  not  included  in  these  earlier  studies.  The 
new  features  are _ad  hoc  methods  in  which  J- shift  and  lasing  cutoff  on  a 
particular  band  are  approximated.  Model  predictions  are  compared  with  those 
of  the  comprehensive  model  of  Ref.  9.  For  cases  of  practical  interest, 
excellent  agreement  is  found.  In  Section  III,  the  usefulness  of  the  model  for 
parametric  variation  is  illustrated.  In  Section  IV,  model  predictions  are 
compared  with  selected  experiments. 
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II.  THEORY 
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A.  MODEL  FORMULATION 

A simple  model  for  a pulsed  HF  laser  pumped  by  the  + E^  chain 
reaction  is  presented.  Initially,  the  medium  is  assumed  to  be  a homogeneous 
mixture  of  F^,  F,  H,  and  diluent.  Relative  concentrations  of  F/F^  and 

H/H9  are  chosen  to  simulate  initiation  by  an  external  source.  The  dominant 

2 10 
reactions  of  the  chemical  system  were  previously  identified  as: 

1.  The  + F^  chain: 

gc(v)kc 

F + H2 ► HF(v)  + H 

8h(Y)kh 

H + F2 ► HF  (v)  + F 

2.  Vibration  - translation  (VT)  deactivation: 

khf(M ’ 

HF(v)  + M HF(v-  1)  + M 

3.  Vibration  - vibration  (VV)  transfer: 


HF(v)  + HF(v') ►HF(v-l)  + HF(v'+l) 

k / 
vv 

HF(v)  + H2(v7) ► HF(v-l)  + H2(v'  + 1) 


Dissociation  - recombination: 


kdf(M) 


f2  + m- 


F + F + M 


kdh(M) 


Hz  + M 


H + H + M 


kdhf(M) 


HF(v)  + M- 


•H+FtM 


The  overall  rates  of  the  chain  reactions  h and  lc,  and  the  fractions  of  the 

c h 

population  produced  in  the  reaction  product  HF(v)  given  by  g^(v)  and  g^(v)  are 
easily  obtained  from  the  rate  coefficients  of  Reactions  (1)  and  (2)  in  Table  I. 
For  the  constant-gain  method  of  solution  used  in  this  model, 
it  is  lcnown  that  the  effect  of  the  (VV)  reactions  in  item  3 above  can 
be  neglected.  Initial  reactant  mixtures  with  total  pressure  less  than  2 atm 
are  considered;  therefore,  the  di ssociation- recombination  reactions  in  item 
4 are  neglected. 

The  rate  of  change  of  HF(v)  may  be  written  as^ 


d[HF(v) 

dt 


= P+  X-  X , + D 

V V V-  1 


1 - ,u 

V+  1 , V V,  V-  1 


where  chemical  pumping  is 


pv  = gc(v>kc[F]  rn2]  + gh(v)kh[H]  [f21 


vibrational  deactivation  is 


Dv,v-1  = [HF(v>]  £khf  (M’V)  fM] 


and  Xv  is  the  rate  of  laser  photon  emission  for  the  band  with  lower  level  v. 
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Table  I. 


Rate  Coefficients  for  H 2 


Chemical  Laser 


a 


Reaction 

No. 

Reaction 

Rate  ( ocffiCient 
icrn^fmolc-src ) 

M. 

h 

la 

F ♦ H, 

10)  II F ( 11  ♦ II 

kL 

Z 6 > 

,o"  ' <■« 

lh 

F ♦ H, 

(0)  MF(2)  ♦ H 

klb 

tc 

F * 

10)  HFI  3>  4 H 

k|c 

4.  4 « 

,o"-‘ 

Id 

F ♦ H, 

(0)  HF(4>  4 H 

k Id 

le 

F 4 H- 

(0)  HFI  5)  ♦ H 

k-U- 

1.1 

If 

F ♦ H2 

(0)  H F 1 6 1 4 II 

k-l( 

1.9. 

io' »-°  S'  s 

2*» 

H*  L 

HF(0)  ♦ F 

ki. 

Zb 

H 4 F2 

HFI  1)  » F 

Lb 

2.  S . 

,0 

2c 

H ♦ F, 

HFI 21  4 F 

k , 

2d 

H * F2 

HFI  3)  ♦ F 

k2«f 

2e 

" * f'z 

HF(4)  ♦ F 

k , 

Zf 

H 4 F2 

HFI  S)  4 y 

Lf 

2c 

H * F, 

HFI6I  4 F 

k , 

4 8 . 

,0I*.2.4.5 

2h 

H 4 F, 

HFI  7)  4 F 

Lh 

S.  5 X 

10  1 * *"  4 6 

2i 

H ♦ Fj 

HFI  8)  4 F 

k2, 

2.  S « 

|012-2  4 5 

JavC 

HF|v) 

4 Mj  HFlv.l)  4 M j 

k<- 

v(IO 

* ■ 10°  L 

'■4, 

V,l 

HF,  v 

1 8 

HF«I) 

> HFI0)  4 M, 

k W>| 

1 . *i  ■ 

8 10  ' T 

M2 

3b, 

HFI2) 

4 M HFIII  ♦ M, 

Lb, 

3b 

HF|vl 

♦ M,  HFIv-ll  4 M2 

k<b 

1.5 

. I0I0T 

v 

3 . . 8 

3c 

HF(v> 

4 M , HFIv  - 1)4  M . 
4 4 

kuv 

18  * 

I0*4T4)v 

M4 

= F,;  v 

1 . .8 

3d 

HF|v| 

* Ms  HFlv.l)  » 

kb.y 

v(8. 

7 „ 10’  't'| 

Ms 

lie;  % . 

1 . 8 

3c 

HF|v) 

♦ M,  HFIv'l  4 Mh 

k.lcv 

1.8 

. io1'"0  7/9 

M 

H.v  1 

8.  v’  • v 

3(v 

HF(vl 

4 M?  HFIv-  1)  4 M? 

k,fv 

v(  1 * 

iovr2) 

M, 

Hy  v 

1 8 

4a 

HF|v| 

4 HFIv)  HFIv-  1)  4 HF!v4  1) 

k 4- 

1 . 5 * 

10 1 “T* ' “ 

v 

1 . 7 

4b 

HF|v| 

4 HF|v4  1 1 HFIv-')  4 HF(v42) 

k4h 

0.  Sk  . 

4a 

» 

1 . b 

4c 

HF|v) 

4 HF(v42)  HFIv-  II  4 HF(vf  3) 

L. 

0. 2Sk 

4a 

V 

1 . .5 

4d 

HF(v) 

4 HFIv*  3)  HFIv-  1)  4 HF(v44l 

k4d 

0.  USk 

1 . . 4 

Sa 

HF(Q| 

♦ h2id  HFII)  * H2(0) 

L. 

9 . 10  1 * 

Sb 

HF(  1) 

4 11,111  H F 1 2 ) 4 H2<0) 

Lb 

2.9* 

io12 

Sc 

HF(2) 

4 ||,in  HFI  31  4 H2<0| 

Lc 

9 x 10 

Sd 

HF(  3) 

4 H2<ll  HFI4)  » H2(01 

Ld 

2 • I0n 

Sc 

HF(O) 

4 H2I2)  HFII)  4 H2(l) 

L.. 

L. 

Sf 

HFCII 

4 11,(2)  HFI2I  • H2II) 

Lf 

Lb 

f>«v 

M2<v» 

4 M8  ll2(v  - 11  4 Mm 

k 

f'A  v 

v|2 . 

■.  . IO'4, T4  ' 

' 2'M* 

all  eweept  H.  H, 

7b 

H2(v|  < 

i M„  II/'  II  I M„ 

k6bv 

v(  10 

V 'l 

v 

H.  H£ 

Not*-  that  dt*s»«  iatn  n r«*<  onihmai  > 1 reactions  hav.  ’>■  . n m lr«  ted 


5 ■ 

M r<illi«i'in  partnrr,  v vibratxmal  l«*v«>|. 


If  we  assume  that  the  highest  level  pumped  by  the  chain  reaction  is  HF(N), 

and  that  lasing  occurs  on  HF(N) ►HF(N-  II  and  all  lower  bands,  Eqs.  (1) 

may  be  solved  for  to  yield 


The  total  rate  of  increase  in  the  HF(v)  population  is  related  to  the  pumping 
rate  by 


N N 

V)  d[HF(v|]  p 


The  integral  of  this  relation  may  be  stated  as  (at  time  t) 


E HF(''»  = 2|lH21e=o-lH2'1| 

v -0 


(3a) 


(3b) 


The  total  rate  of  photon  emission  is  given  by 


N-  1 


The  refore , 


X 


T 


N-  1 

S °V+1’V 


(4) 


(3) 


During  lasing,  we  assume  that  the  gain  at  line  center  for  the  transition 
lasing  in  each  band  is  maintained  at  the  cavity  threshold  value  a ^ 


Ot 


thr 


1 

2 L 


ln(RQR 


(6) 


where  L is  the  length  of  the  active  medium,  Rq  is  the  reflectivity  ol  the 
output  coupler,  and  R accounts  for  a combination  of  the  reflectivity  of  the 

L-j 

second  mirror  and  other  cavity  losses.  The  gain  at  line  center  on  the  P-branch 
of  the  transition  with  lower  level  J is 


a ( v , J ) 


a)  (v,  J)  0 (v,  J)  j 


f(2  J+  1 ) 
2 J - 1 


HF(v+  1 , J-  1)  - H F ( v , J ) 


(7) 


where  is  Avogadro's  number,  h is  Planck's  constant,  >v(v,  J)  is  the 

wavenumber  of  the  transition,  o(v,  J)  is  the  Voigt  line  shape  function  evaluated 

at  line  center  (Appendix  B),  and  the  Einstein  isotropic  absorption  coefficient 

is  B , = 6.2828  X 10^°  M , . | 2 J/(2J+  1),  cm2 /molecule  -J-sec.8  We 

vj  VT 1 , V g 

approximate  the  electric  dipole  matrix  element  by 


M , = (v+  1 ) 1 /Z  x [0.988  x 10"19  X (1  + 2.66  x 10~2J)] 

v+  I , v 


The  value  of  o(v,  J)  is  taken  as 


8 


(8) 


to  (v,  J) 


60.  32  T 

^v  ~ TTTWJ 


(9) 


where  the  o used  are  given  in  Table  II. 
v 

A Boltzmann  distribution  of  the  rotational  levels  is  assumed,  and  trans- 
lational and  rotational  degrees  of  freedom  are  taken  to  be  in  equilibrium;  therefore 

e (v) 

HF(v,  J)  = HF(  v)  — -7p-  (2J+  1 ) exp[  - J (J  + 1 ) 9r  (v)  / T ] , (10) 
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Table  II.  Band  Energy  Level  Spacings 


T ransition 
(v+1),  V 

u)  (cm  ) 

V 

1.0 

3950 

2.  1 

3790 

3.2 

3630 

4.  3 

3470 

5.4 

3320 

6 , 5 

3170 

7,6 

3020 

8.7 

2880 

o 

OD 

27  30 

-14- 


where  9 (v)  is  approximated  by 

9 (v)  ( 30  - v)K  . (11) 

r 

Because  o(v,  J)  a ^ during  lasing  on  the  v+ 1 *-v  band,  Eq.  (10) 

may  be  substituted  into  (7)  to  determine  a relationship  between  populations  in 
the  levels  of  lasing  bands.  Therefore,  dviring  lasing,  the  concentrations 
satisfy  the  relation 


HF(v+  1 ) - z(v,J)  + HF (v)  T1 7 (v) 


(12a) 


whe  re 


z ( v , J ) 


4'trTo' 


thr 


hN^o(v7J)  0(v,  J) 


J(J-  1) 

Y v+1 

b : 9 (v+  n (2.i+  n 

vj  r 


(12b) 


and 


Vv> 


9r(v) 

TTv+T) 

r 


J(J-l) 

Yv+1 

J(J+1) 

Yv 


(12c) 


and 


9 (v)/T 

V = e r . ( 1 2d) 

1 v 

The  time  history  of  the  temperature  may  be  determined  from  the  energy 
equation 
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dT  _ 

. cTT  = 
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V'd|N.1 

■pl-E^tt 
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(13) 
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where  H is  the  molar  enthalpy,  C is  the  molar  specific  heat  at  constant 
> i 

volume  of  species  i,  and  P is  the  laser  power  per  unit  volume  defined  by 


N-  1 

PL  = hcNA]C  u)(v,J) 


(14) 


v=0 


B.  APPROXIMATE  MODEL  SOLUTION 


An  approximate  method  for  solution  of  the  model  equations  is  discussed. 
This  method  incorporates  various  physical  approximations  to  facilitate 
solution  by  simple  numerical  computation.  Recall  that  we  have  already 
included  the  following  approximations: 

The  mixture  is  assumed  homogeneous,  and  initiation  is 
simulated  by  introducing  an  initial  F and  H atom  concentration. 

2.  During  lasing,  the  gain  at  line  center  is  maintained  at  the 
threshold  value. 

3.  Vibration-vibration  relaxation  mechanisms  may  be  neglected. 

4.  The  rotational  levels  are  in  equilibrium. 

5.  Reverse  reactions  may  be  neglected. 

In  order  to  obtain  a simple,  but  useful,  solution,  some  additional 
assumptions  and  approximations  are  made.  With  operation  of  the  chain,  the 
rate  of  change  of  the  F-atom  concentration  is  given  by 


[H]  [F?l  - [F]  [H21  k 
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and  the  simple  relation  between  the  chain  carriers  gives 


[F  + H]t  = [F  + H]t; 


(15) 


(16) 
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From  earlier  work,  10  it  is  known  that  the  ratio  [F]/[H]  quickly  assumes  a 
value  given  by  the  steady-state  chain  condition 

[F]  [H2]  kc  ~ [H]  [Fz]  kh  . (17) 

Because  this  ratio  is  approached  on  a time  scale  that  is  short  in  comparison 
with  the  laser  pulse  length,  we  assume  that  Eq.  (17)  holds  for  all  time.  These 
approximations  permit  simplification  of  Eq.  (lb)  to 


Pv  = Gvkc  (F]  [H2] 

whe  re 


(18a) 


Gv  = gc(v)  + gh(v) 


(18b) 


By  approximating  the  heat  released  to  be  that  of  the  chain,  Eq.  (13)  becomes 


£ fNi1cv.  w ~ IP]  Ih2]  kc AH  - PL  - h,E''  • 


(19) 


where  AH  is  129 . 6 kcal /mole , and  hf  is  taken  to  be  1 1 . 0 kcal /mole. 

The  number  of  monatomic  and  diatomic  molecules  is  not  changed  by  the 
reaction;  therefore,  we  assume  that  the  specific  heat  at  constant  volume  of  the 
mixture  is  constant  and  equal  to  its  initial  value. 

To  facilitate  determination  of  approximate  HF(v)  concentrations,  we 
further  assume  that  (1)  the  lasing  bands  terminate  in  sequence  with  the  1-0 
band  shutting  off  first  and  that  (2)  all  bands  are  lasing  on  transitions  that 
have  the  same  value  of  J.  Termination  is  assumed  to  occur  when  X , as 
given  in  Eq.  (2),  first  becomes  negative.  Previous  studies  with  the  com- 
prehensive model  of  Emanuel  indicated  that  the  1-0  band  shuts  off  first  in  the 
pulse,  whereas  the  other  bands  shut  off  nearly  simultaneously  at  some  later 
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time,  or  those  that  lase  for  longer  times  do  so  at  greatly  reduced  intensities, 
in  Hence,  essentially  all  higher  bands  terminate  at  nearly  the  same  time  and  at 

a time  somewhat  later  than  for  the  1-0  band.  This  behavior  is  accurately 
predicted  by  assumption  (1)  and  our  somewhat  ad  hoc  definition  of  termination. 

One  cannot  expect  a gain-equals-loss  model  to  predict  time- resolved 

14 

spectra  such  as  given  by  Hough  and  Kerber.  Therefore,  the  second 
assumption  (2)  that  all  transitions  lase  on  the  same  value  of  J is  not  a severe 
one.  However,  the  fact  that  this  J must  shift  during  the  run  in  order  to  predict 
accurate  results  is  significant.  We  assume  a value  of  J that  maximizes  the 
population  of  the  lower  level  of  the  lowest  lasing  band.  This  implies  that  the 
excited  HF  population  approaches  the  minimum  necessary  to  maintain  lasing 
on  these  bands.  This  result  is  linked  to  the  constant-gain  assumption  and 
is,  in  fact,  compatible  with  it.  Typically,  the  J values  that  occur  are  be- 
tween 2 and  17.  This  J- shift  is  necessary  for  the  calculations  to  yield  time 
histories  of  HF(v)  that  are  comparable  with  those  of  the  comprehensive  model. 

These  assumptions  lead  to  some  notation  simplifications;  for  example, 

we  let 


z = z(v, J) 
v 


(20) 


and  similar  definitions  for  y , x , 0 , and  B are  made.  Therefore, 

‘v  v v v 

Eq.  (12b)  takes  the  form 


z 

v 


, _ J ( J - 1 ) 

4tt  Tor,,  y , , 
thr  v + 1 

hN  . u)  0 B S (v+  1 ) (2  J+  1 ) 
A v v v r 


(21) 


and  Eq.  (11a),  which  relates  populations  associated  with  lasing  bands,  becomes 


HF  (v+1)  z + HF(v)T) 
v V 


(22) 
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li  wo  denote  by  v-'-  the  lower  level  of  the  lowest  band  lasing  at  a given  time, 
i'  can  write  the  rate  of  change  of  the  nonlasing  HF(v)  populations  as 


v -0 


v - 1 

H F ( v ) V 
v-Q 


P + D . 

V v:  , 


V : - 1 


(23) 


From  Eq.  (22),  which  holds  for  v > v"  , we  can  show  that 


N 
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X.  II  T 

V— V v ,=v/  + 1 

The  integral  of  Eq . (23),  when  summed  with  Eq.  (24),  gives  the  sum  of  HF(v) 
over  all  v,  which  is  also  given  by  Eq,  (3b). 

Therefore,  we  can  determine  HF(v  ) as  given  in  Eq.  (24)  by  combining 
these  relations.  The  preceding  equations  can  be  solved  numerically  for  the 
time  histories  of  HF(v),  T,  Xv<  etc.  provided  that  a value  for  J is  specified. 
(We  use  a simple  forward  difference  routine.)  At  each  integration  step,  we 
choose  the  value  of  J that  makes  HF(v't)  a maximum.  This  population  is 
normally  driven  to  a maximum  as  a result  of  stimulated  emission  from 
radiation  fields . 


COMPARISON  WITH  COMPREHENSIVE  MODE  I. 


In  Fig.  1,  time  histories  of  intensity  for  various  levels  of  initiation 
are  compared  with  the  predictions  of  the  comprehensive  model  of  Emanuel 


et  al.  The  figure  shows  that  our  model  does  quite  well  for  [F  + H]/|F?] 

0.  1 and  0.01.  However,  for  f F + Hj/jF^]  0.001,  we  find  that  the  simple 
model  is  somewhat  low,  and  that  the  comprehensive  model  gives  a surprising 
result.  The  spike  near  the  end  of  the  calculation  with  the  comprehensive 
model  was  caused  by  production  of  F atoms  from  thermal  dissociation  owing 
to  the  increase  in  temperature.  This  phenomenon,  which  is  dependent  on  the 
bond  strength  of  F7,  should  be  viewed  with  some  skepticism  until  experimental 
corroboration  is  obtained.  Since  the  present  model  does  not  include  thermal 
dissociation,  a second  calculation  with  the  comprehensive  model  was  made 
that  neglected  F^  dissociation  [dotted  curve  of  Fig.  1(c)].  With  dissociation 
suppressed,  the  agreement  between  the  models  is  greatly  improved.  As 
stated  earlier,  the  primary  goal  of  this  work  is  to  develop  an  efficient 
method  of  predicting  pulse  energy  and  time  histories  of  pulse  intensity.  We 
find,  from  Fig.  1,  that  the  predicted  intensities  are  closely  approximated 
except  near  pulse  termination.  At  termination,  the  simple  model  predicts  a 
more  abrupt  shutoff  than  the  more  extensive  calculations  indicate.  However, 
the  cutoff  time  of  the  simple  model  is  an  excellent  indication  of  the  time  to 
release  most  of  the  pulse  energy.  Pulse  energies  for  the  two  models  are 
compared  in  Table  III.  Over  the  wide  range  of  initiation  levels  considered, 
the  pulse  energies  from  the  two  models  agree  to  within  5%. 

The  step-size  sensitivity  of  the  simple  model  was  examined  by  varying 
the  number  of  steps  used  in  determining  the  time  history  of  HF(v)  between 
ino  and  1000.  The  answers  varied  by  less  than  3%  for  most  cases  studied. 

The  computational  time  of  the  simple  model  is  approximately  1 . 0 to  1 . 3%  of 
the  comprehensive  model.  In  fact,  all  200  simple-model  calculations 
carried  out  in  this  study  were  performed  with  less  computer  time  than  was 
required  for  two  typical  cases  with  the  comprehensive  model. 
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Tabic  III.  Comparison  of  Pulse  Energy  as  Predicted  by 
the  Simple  and  the  Comprehensive  Model 


Comprehensive  Model 

Simple  Model 

[f+h]/[f21 

Pulse  Energy 

( J /liter ) 

0.  1 

291 

306 

0.0  1 

208 

205 

0.00  1 

103  ( 1 0 0 ) a 

100 

Calculation  with  dissociation  of  F?  suppressed. 


Hi.  PARAMETRIC  STUDY  OF  LASER  PERFORMANCE 


In  this  section,  the  effect  of  varying  the  initial  gas  mixture  and  cavity 

conditions  is  examined.  The  laser  pulse  is  characterized  by  its  duration  t , 

peak  power  P , and  total  energy  E.  The  variations  are  made  about  a standard 

cast.-  with  1 F-,:1  II,:0.01  (II+F):20  He  with  T equal  to  300  K and  P.  equal  to 

1 atm.  The  standard  cavity  conditions  have  R , R.  , and  L equal  to  0.9,  1.0, 

and  100  cm,  respectively.  For  brevity,  only  in  Fig.  2 do  we  include  results 

for  peak  power  P and  the  time  at  which  it  occurs  t . 

P P 


A. 


THRESHOLD 


As  expected,  we  find  that  E and  t increase  with  decreasing  cavity 

losses.  These  trends  are  shown  in  Fig.  2.  The  agreement  between  the  two 

models  on  the  prediction  of  pulse  energy  is  very  good  for  the  case  [F+H]/ 

[F  ] 0.01.  Also  shown  in  Fig.  2(a)  is  a family  of  curves  obtained  from  the 

simple  model  for  several  different  levels  of  initiation.  The  predictions  of 

P also  compare  favorably  with  the  comprehensive  model  shown  in  Fig.  2(c). 

As  shown  in  Fig.  2(b),  the  two  models  agree  very  well  on  predictions  of  t ; 

however,  the  agreement  of  values  of  t are  poorer.  (This  result  is  also 

shown  in  Fig.  1.)  This  is  because  the  comprehensive  model  tends  to  yield 

pulses  with  long  low  power  tails,  whereas,  the  simple  model  shows  a more 

abrupt  cutoff.  Except  for  this  discrepancy,  the  models  are  found  to  agree 

very  well.  The  behavior  observed  here  is  similar  to  that  observed  earlier  for 

the  effect  of  cavity  loss  on  ail  pulse  characteristics.  As  expected,  the 

laser  performance  is  more  sensitive  to  cavity  threshold  for  low  levels  of 

initiation  [Fig.  2(a)],  This  may  be  noted  by  observing  the  change  in  the 

-4 

character  of  the  curve  for  X = 10 
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B . 


INITIATION  LEVEL 


The  initial  atom  concentration  was  varied  while  all  other  conditions 

were  held  fixed.  Again,  the  behavior  predicted  is  very  similar  to  that  found 

earlier,  ^ except  that  we  are  able  to  examine  efficiently  the  behavior  at  very 

low  levels  of  initiation  with  the  simple  model.  The  results  are  shown  in 

Fig.  3.  The  decrease  in  t and  the  change  in  slope  in  the  curve  for  E as 

C -4 

[ F+I I] /[F  , ] decreases  below  10  reflect  the  sensitivity  of  long-chain  systems 
to  cavity  losses  (threshold).  The  agreement  between  the  comprehensive  model 
and  the  simple  model  in  predicting  E over  a wide  range  is  very  good  [Fig.  3(a)]. 
The  times  used  for  the  comprehensive  model  in  Fig.  3(b)  are  the  times  at  which 
the  power  has  dropped  to  10%  of  its  peak  value.  Examination  of  Fig.  1 shows 
that  this  is  a better  indicator  of  pulse  length  than  is  t . It  was  desired  to 
compare  with  the  comprehensive  model  at  even  lower  values  of  [F+H]/[F-,]; 
however,  the  computations  with  it  became  prohibitively  long  in  that  regime. 


CL  DILUENTS 

In  Fig.  4,  the  effect  of  varying  the  helium  diluent  concentration  is 
examined  for  several  initial  mixture  pressures.  Here,  as  the  diluent  is 
increased,  reactants  decrease  because  total  pressure  is  held  constant.  The 
results  are  consistent  with  our  expectation  that  increased  diluent  at  the  ex- 
pense of  reactants  would  yield  lower  energies.  Also,  longer  pulse  lengths 
result  because  of  lower  reactant  pressures  and  temperature.  The  dashed 
curves  illustrate  the  effect  of  diluent  on  pulse  characteristics  for  mixtures  of 
constant  reactant  concentration.  Since  the  VT  deactivation  rate  by  He  is 
small  compared  with  that  of  HF  by  itself,  this  result  is  primarily  an  effect 
of  temperature  control.  Again,  the  agreement  is  excellent  between  the  two 
models  on  the  predicted  pulse  energy. 
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SIMPLE  MODEL 

COMPREHENSIVE  MODEL 

SIMPLE  MODEL  CURVE  OF  CONSTANT 
REACTANT  CONCENTRATION 


Effect  of  He  Concentration  for  Selected  Values  of  Initial  Pressure. 
The  dashed  curves  represent  loci  of  constant  reactant  concentration 
Curves  w,  x.  y,  and  z represent  F 2 concentrations  of  505,  252,  126, 
and  54  Torr,  respectively,  a.  Effect  on  pulse  energy,  b.  Effect  on 
pulse  duration.  Initial  gas  mixture:  F + H : H2  : F2  : He,  0.01  : 1 
: 1 : He;  Tj  300  K;  cavity  conditions:  Rq  0.9,  1-0,  L 

100  cm. 


D.  INITIAL  CONCENTRATION 

The  initial  [P^]  was  varied  while  other  conditions  were  the  same  as 
those  for  the  standard  case  for  typical  levels  of  initiation.  The  pulse  energy 
shows  a maximum,  and  the  pulse  duration  shows  a minimum  when  the  mixture 
is  nearly  stoichiometric  (Fig.  5).  The  behavior  as  [H^]  increases  to  the 
stoichiometric  mixture  is  clear;  more  reactant  increases  the  reaction  speed 
and  energy  released.  For  mixtures  significantly  above  stoichiometric,  we 
find  that  [H^]  is  increasing  at  the  expense  of  [F^];  therefore,  we  observe  a 
decrease  in  pulse  energy  and  an  increase  in  pulse  duration.  The  curves  are 
fairly  flat  in  the  region  just  above  stoichiometric,  which  indicates  a very 
weak  dependence  of  pulse  characteristics  on  [H^]  in  this  region.  This  be- 
havior is  observed  over  a wide  range  of  initiation  levels  (Fig.  3).  Again, 
the  simple  model  yields  pulse  energies  in  excellent  agreement  with  the 
comprehensive  model. 


E.  INITIAL  TEMPERATURE 


The  initial  gas  mixture  concentration  is  inversely  proportional  to  the 
initial  temperature  for  a constant  pressure  system.  Therefore,  the  behavior 
seen  in  Fig.  6 is  caused  mainly  by  the  amount  of  reactants  available  to 
generate  laser  energy.  Recall  that  the ‘rate  coefficients  for  all  reactions  are 
less  well  known  below  300  K,  and  predictions  in  this  regime  may  be  suspect. 

The  disagreement  between  the  two  models  observed  in  Fig.  6(a)  is  the 
greatest  found  in  the  entire  study.  Above  400  K,  the  two  models  begin  to 
differ  considerably.  However,  this  does  not  appear  to  be  a severe  problem 
as  it  is  not  anticipated  that  any  useful  system  will  be  initiated  at  elevated 
temperatures. 


:V.  COMPARISON  WITH  EXPERIMENT 


Since  this  model  has  been  found  to  be  equivalent  to  the  comprehensive 
model  over  the  wide  region  of  interest  for  most  practical  lasers,  its  ability 
to  predict  energy  and  power  of  experimental  devices  is  comparable  with 
that  of  the  comprehensive  model. 

We  have  used  the  model  in  an  effort  to  predict  the  experimental  results 
obtained  by  Hofland  et  al.  ^ First,  the  model  was  compared  with  experiment 
by  using  the  rate  coefficients  listed  in  Table  I.  The  value  of  [F  + H]/[F9] 
was  selected  by  matching  pulse  energies  at  a selected  value  of  H ^ concentra- 
tion. The  estimated  [F  + H]/[F-,]  ratio  was  0.005  for  cases  with  argon  and 
0.003  for  cases  without  argon.  Although  the  predicted  energy  followed  that 
of  experiment  quite  well  for  variations  in  [ 1 ] , the  predicted  pulse  duration 
was  about  six  times  longer  than  observed  in  the  experiment.  The  value  of 
[F  + Hj/fF^]  was  also  estimated  to  be  0.06  by  matching  pulse  durations.  The 
calculated  pulse  energies  were  then  roughly  a factor  of  4 larger  than  observed 
experimentally. 

Finally,  the  rate  coefficient  for  HF-HF  V-T  relaxation  was  assumed 
to  be  of  the  form 


k 


3 

av 


1 


as  suggested  by  the  recent  results  of  Kwok.  In  Fig.  7,  calculations  using 
our  model  and  this  rate  are  compared  with  experiment.  The  conditions  of 
the  experiments  arc  listed  in  the  figures.  An  additional  loss  of  10%  was 
assumed  to  account  for  contingencies  such  as  volume  and  surface  scattering 
losses.  The  value  for  [F  + Hj/fF^]  was  determined  by  matching  the  pulse 
energies  at  - 3.  For  the  mixtures  without  argon,  [F  + Hj/fF^]  was  0.005, 
and  for  mixtures  with  argon  [F  + Hj/fF^]  was  0.009  and  0.006  for  cases  with 
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[He]  = 36  and  54,  respectively.  On  doing  this,  we  found  that  the  model 
predicts  pulse  durations  within  10%  of  those  observed  experimentally. 

'V  Although  we  expected  discrepancies  between  experiment  and  theory,  we 

emphasize  them  to  illustrate  the  need  for  the  simple  model.  The  simple 
model  may  be  used  ( 1 ) to  test  efficiently  the  sensitivity  to  rate  coefficient 
variations,  and  (2)  to  incorporate  efficiently  the  effects  of  plasma  kinetics 
or  mode  structure.  If,  after  such  modifications,  this  model  begins  to  com- 
pare favorably  with  experiment,  we  can  then  return  justifiably  to  the  more 
expensive  comprehensive  calculations. 


DENSITY  (J/liter) 


E (kV/cm)  Cnn  (%) 


= 1200 


A A 

£ a 


■ * THEORY 

□,o,  a experiment 


PERCENT  H2,  X 

Figure  7.  Comparison  of  Theory  and  Experimental  Results  of 
Hofland  ct  al.6  at  800  Torr.  Note  that  the  calculated 
and  observed  pulse  duration  is  3 psec  for  the  6 F2:3 
H2:54  He:37  Ar  mixture.  For  cases  with  Ar,  the 
value  of  (F  + H)/(F2)  used  in  the  model  is  0.009  and 
0.006  for  cases  with  [He]  = 36  and  54,  respectively, 
and  0.005  for  the  case  without  Ar.  CQp%  is  the 
percent  transmitted  by  the  output  coupler. 

V.  CONCLUSIONS 


"1 


> 

A simplified  model  of  the  + F^  pulsed  chemical  laser  has  been 
presented.  This  model  yields  results  in  agreement  with  a more  compre- 
hensive model  at  about  1%  of  the  cost.  Therefore,  we  expect  this  model  to 
be  the  logical  tool  for  use  in  improving  the  agreement  between  theory  and 
experiment.  Once  this  is  done,  we  can  then  return  to  the  more  detailed 
calculations  of  the  comprehensive  model  (RESALE)  or  of  the  model  given 
in  Ref.  14. 

In  addition,  the  economy  of  making  accurate  theoretical  predictions  will 
permit  the  model  to  be  used  extensively  in  parametric  studies  or  as  a 
"building-block"  for  incorporation  into  other  models  for  examining  effects  of 
mode  structure  or  other  physical  phenomena,  such  as  plasma  kinetics,  on 
the  performance  of  laser  systems. 


APPENDIX  A 


The  detailed  chemical  kinetic  model  used  for  this  study  is  presented 

here.  The  rate  coefficients  given  in  Table  I are  based  on  the  review  of  the 

22 

rate  measurements  and  predictions  given  by  Cohen.  The  entire  kinetic 
model  was  used  for  the  calculations  with  the  comprehensive  model  (RESALE) 
presented  in  this  report.  Only  the  appropriate  reactions  and  rate  coefficients 
were  used  in  the  simple  model. 
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APPENDIX  B 
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The  Voigt  profile  at  line  center  is  given  by 


2 

ey  [l  -erf(y)]/aDp(v,  J) 


where 


y 


(ln2)  aLR/aDP(v’ J) 


The  Doppler  half-half  width  is 


( v > J ) 


2 X 10  NAkT  ln2 
A 

w'" 


,1/2 


®(v, J) 


(Bl) 


( B2) 


( B3) 


and  the  Lorentz  half-half  width  is 

*LR  = X)  fNi]  RTl/2  W2  aj  [Njl/P  > (B4) 

i j 

where  p is  the  density,  W is  the  molecular  weight  of  HF,  and  W is  the 
molecular  weight  of  the  total  mixture.  Only  billiard-ball  collisions  have  been 
included  in  Eq.  (B4). 
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THE  IVAN  A.  GETTING  LABORATORIES 


r 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility ard  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 

fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 

pheric optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 

propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  appl  ed  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment,  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields,  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  spare  on  space  systems. 

THE  AEROSPACE  CORPORATION 
El  Segundo,  California 


